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Abstract. 

Based on the 1960 idea of Lipkin, the minimization of energy of a symmetry- 
restored mean-field state is equivalent to the minimization of a corrected energy of a 
symmetry-broken state with the Peierls-Yoccoz mass. It is interesting to note that 
the "unphysical" Peierls-Yoccoz mass, and not the true mass, appears in the Lipkin 
projected energy. The Peierls-Yoccoz mass can be easily calculated from the energy 
and overlap kernels, which allows for a systematic, albeit approximate, restoration of 
translational symmetry within the energy-density formalism. Analogous methods can 
also be implemented for all other broken symmetries. 



PACS numbers: 21.10.-k, 21.60.-n, 21.60.Fw 
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1. Introduction 

One of the biggest successes of the mean- field (MF) and energy-density-functional (EDF) 
methods applied to many-body systems consists in including correlations through the 
mechanism of the symmetry breaking. A single broken-symmetry MF state or a one- 
body broken-symmetry density can, in fact, represent a large class of correlations 
that are physically important. Moreover, the breaking of symmetry on the MF level 
corresponds to the appearance or disappearance of certain correlations that in finite 
systems may have all features of the phase transitions. 

A link between the broken-symmetry MF states and correlated symmetry- 
conserving states is provided by the symmetry-restoration methods [1]. The broken- 
symmetry states can be viewed as auxiliary objects, which facilitate obtaining the real 
quantum mechanical states that have all good quantum numbers. Within the EDF 
methods, this leads to a very fruitful idea of the projected energy being a functional 
of the symmetry-breaking one-body density [2, 3]. However, variational equations then 
become quite difficult to solve, and the full symmetry-restored MF or EDF methods 
were to date applied only in several particular cases [2, 3, 4]. 

Numerous approximate methods to restore broken symmetries by variation before 
or after projection (VAP) were formulated and implemented in the past [1, 5]. The 
main issue here is the feasibility of the method - it must be a reasonable compromise 
between the precision and closeness to the exact VAP on one side, and with the numerical 
effort to execute it on the other side. In this respect, the Lipkin method [6] has very 
many advantages, which are discussed in the present study. It allows for calculating an 
approximate VAP energy without any necessity to perform the projection at all. In the 
past, it has been mostly used for the particle-number restoration, within the so-called 
Lipkin-Nogami [7] formulation. However, the original Lipkin method was formulated 
for restoring the translational symmetry, and in this study this case is studied in detail. 

2. Results 

To fix the notation, let us begin by recalling the basic standard definitions and properties 
pertaining to the translation symmetry. Let |$) denote a normalized Slater determinant 
built of single-particle orbitals that are localized in space. Since the total momentum 
operator P = J2iLi Pi is the generator of translation, states |$) can be shifted by R to 
an arbitrary location in space as 





(2) 
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and P\P) = P\P)- The normahzation condition of states \P) is chosen in such a way 
that the original Slater determinant is a simple integral thereof, 

1$) = |$(o)) = j d?P\P). (3) 

The Slater determinant |$) is, therefore, a normahzed wave packet built of 
non-normalizable center-of-mass plane waves |P). For a system described by a 
translationally invariant Hamiltonian [H, P] — 0, one can determine the average 
energy E^{P) of each plane wave, which is called the projected energy, as 

E^{P) = ($|^|P)/($|P). (4) 



2.1. Kernels 

Prom Eq. (2) one sees that the projected energy (4) is given by the Fourier transforms 
of the overlap and energy kernels, which are defined by 

I{R) =($|$(P)>, (5) 
H{R)^{^H\^R)), (6) 

respectively. Therefore, properties of the kernels must be discussed first. Moreover, all 
results below depend only on the kernels; hence these results apply automatically to the 
EDF approaches, where very often does not start with the Hamiltonian but the diagonal 
energy density is extended to the energy kernel [8]. 

In Figs. 1 and 2 are shown, respectively, logarithms of the overlap kernels, ln(/(i?)), 
and reduced energy kernels, h{R) = H{R)/I{R), calculated in 9 doubly-magic spherical 
nuclei from "^He to ^°^Pb. Calculations were performed for the SLy4 [9] parametrization 
of the EDF, by using the code HFODD [10] (v2.40g) and the harmonic-oscillator (HO) 
basis of up to A^o = 14 shells. 

It can be seen that in the case of the translational symmetry, the so-called Gaussian 
Overlap Approximation (GOA) [1], given by 

I{R) = exp(-|ai2'), (7) 

h{R) ^ho- \h2R^, (8) 

is excellent. Sudden deviations of h{R) from the parabolic dependence on R, which 
can be seen in Fig. 2 around \R\ = 4fm, are due to the finiteness of the HO basis used 
in the calculations. Nevertheless, values of /i2 can be very precisely determined in the 
parabohc region. This was confirmed by repeating the calculations for A^o = 20 HO 
shells, whereby the above deviations appear around \R\ — 4.5 fm and the values of /i2 
stay exactly the same. 



2.2. Projected energies 

Within the GOA, the Fourier transforms above can be analytically calculated, and the 
projected energy of Eq. (4) reads [1] 

£^f°^P) = /^o-£;zPM + ^, (9) 
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Figure 1. (Color online) Logarithms of the overlap kernels, ln(/(i2)), calculated in 9 
doubly-magic spherical nuclei (dots). Thick lines represent parabolic fits to the results, 
determined for \R\ < 3fm. 



where the so-called zero-point-motion correction and Peierls-Yoccoz (PY) mass are given 
by 

EzPM = 3h2/2a, (10) 
MpY = n^a^/h2. (11) 

For completeness, in Figs. 3 and 4, I show the mass dependence of the GOA 
parameters a, -Ezpm, and Mpy, along with the fits of the power-law dependencies. 
One clearly sees that the PY mass is not equal to the total mass mA of the nucleus. 
However, as discussed below, it is the PY mass, and not the physical mass mA, which 
is important for the translational-symmetry restoration. 

The average energy of the system at rest, i?$(0), depends on the Slater determinant 
1$), and in what follows we are interested in minimizing this energy with respect to |$), 
that is, in performing the YAP or variation after symmetry restoration. To this end, 
in this study I follow the seminal idea by Harry Lipkin [6], who realized that the YAP 
calculations can be very easily performed by fiattening the function Eq,{P). 

Indeed, from Eqs. (3) and (4) one obtains the so-called sum-rule property, 

j <1^PE^{P){^P) = ($1^1$), (12) 



Lipkin translational- symmetry restoration in the mean-field and energy- density- functional methods^ 




2 4 6 

R[fm] 

Figure 2. (Color online) Same as in Fig. 1 but for the reduced energy kernels, 
h{R) = H{R)/I{R). To guide the eye, thin lines connect calculated values (dots). 

which tells us that the MF energy is equal to the average of the projected 

energies E^{P) weighted by amplitudes Therefore, minimization of 

that is, the standard MF method, corresponds to an entangled minimization of both the 
projected energies E^{P) and amplitudes whereas the VAP method pertains to 

the minimization of the energy -^^(O) only, and, of course, disregards the amplitudes 
{^\P) completely. 

Note that the GOA amplitudes ($|-P) are all positive, 

(*|J')=(f?)""exp(-fif'% (13) 

and hence the MF minimization is bound to underestimate the momentum spread of the 
Slater determinant. It is then obvious that unless the center-of-mass and internal degrees 
of freedom are exactly separated, like is the case for the closed-shell HO [11, 12, 13] or for 
the coupled-cluster states [14], the MF and VAP Slater determinants can be different. 

2.3. The Lipkin method 



The Lipkin idea of flattening the function E^{P) is implemented in the following way. 
Guided, by the GOA result in Eq. (9), one defines the Lipkin operator 

k = kP^, (14) 
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Figure 3. (Color online) Dots show parameters a (bottom) and /12 (top) that are 
calculated by the parabolic fits to the logarithms of the overlap kernels (Fig. 1) and 
reduced energy kernels (Fig. 2). Lines give estimates of the power-lavif dependence. 

for which the projected average values K{P) are defined as: 

K{P) = {^\K\P)/{<!>\P), (15) 
which in turn fulfill the sum rule: 

J d^PK{P){^\P) = (<i>|i:|<i>). (16) 
Then, by subtracting Eqs. (12) and (16), one obtains 

J d^P {E^{P) - K{P)) (<I>|P) = {^\H - k\^). (17) 

It is clear that the fiattest difference Eq,{P)—K{P) is obtained by adjusting the constant 
k so as to best fulfill the equation 

K{P) = E^{P)-E^{{i). (18) 

Note that since -ft'(O) = 0, the exact projected energy of the system at rest, -E#(0), 
must, by definition, appear on the right-hand side of Eq. (18). 

In case that Eq,{P) grows exactly parabolically, one obtains that 

E$(0) = ($1/7-^:1$), (19) 
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Figure 4. (Color online) The GOA zero-point-motion corrections i?zPM and PY 
masses Mpy (inset) calculated for the values of parameters a and ft,2 that are shown 
in Fig. 3. 



that is, the projected energy of the system at rest E^{0) can be calculated without 
performing any projection at all. 

The essence of the Lipkin method is in finding a suitable value of the correcting 
parameter k, which for each Slater determinant |$) must describe the parabolic growth 
of the function £'$(P). We note here that this growth has nothing to do with the 
physical translational motion of the system boosted to momentum P, in which case the 
energy must grow as E{P'^) = U^P'^ /2mA. The function Eq,{P) simply characterizes 
the distribution of projected energies within the Slater determinant, that is, the degree 
of the symmetry breaking in the Slater determinant at rest. Therefore, the correcting 
parameter k has no obvious relation with the true translational mass of the system. 
Moreover, the correcting parameter must depend on all kinds of approximations or 
space truncations, which are made when obtaining the Slater determinant |$), that is, 
the Lipkin method corrects for these approximations too. 

The right-hand side of Eq. (19) can be minimized by the standard self-consistent 
method, whereby the optimum state |$) can be found. At each stage of the iterative 
procedure one has to determine A;, that is, k must parametrically depend on |$). Note 
that at any given iteration of the self-consistent method, this parametric dependence 
must not he varied. Finally, after the iteration converges, -E$(0) is given by the obtained 
minimum value of the Lipkin projected energy (19). Obviously, the quality of the result 
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crucially depends on the quality of the calculation of which now will be discussed. 



2.4- The Lipkin operator 

A practical method to calculate the correcting constant k from Eq. (18) must not involve, 
of course, the exact evaluation of the function E^{P). A method to probe the function 
Ei^{P) without evaluating it explicitly can be formulated as follows. One first remarks 
that the projected energy (4) can also be calculated as: 

E^{P) = mHUP)\P)imn{P)\P), (20) 

which gives the sum rules: 

j d?PE^{P)g^{P){^P) = ($|^^„(P)|$), (21) 

where gn{P) arc arbitrary functions of P. By using different functions one can probe 
the unknown function E^[P). 

An obvious choice of cjn = -P^"^, which, in fact, has been used within the Lipkin- 
Nogami method [7] to restore the particle number, requires dealing with impractical 
many-body operators. A much better option is provided by the shift operators 
[cf. Eq. (1) and (2)], 

^„ = exp(fi4-P), (22) 

defined for a suitably selected values of shifts -R„. Then, the average values on the 
right-hand side of Eq. (21) become equal to the energy kernels H{Rn), which are very 
easy to calculate. 

In practice, the method works as follows. Suppose one wants to evaluate the Taylor- 
expansion coefficients of -E$(P) up to a given order of 2M, 

M 

E^P) = J2 4'"^-P'"- (23) 

m=0 

After inserting this into Eq. (21), one obtains the set of linear equations, 

M 

Y.A^^eS'^^ ^H(R^). (24) 

m=0 

where the matrix A„„ is defined by the kernels of the momentum operators: 

Anm = P2m{Rn) = ($|P'"^|$(i^„)) . (25) 

In the simplest case of the quadratic Lipkin operator, Eq. (14), one only needs the 
expansion up the second order, that is, for M = 1. Then, by using two points Rq = 
and Ri = R one obtains 

A=( ' ' ^^l^'l^^ ^ (26) 
\ i^MR)) , ($|p2|$(i2)) )' ^ ' 
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This matrix can be easily inverted and then one obtains the first two Taylor expansion 
coefficients E^^ and which are required in Eq. (18), that is, 

h{Q)p2{R) - h{R)p,{Q) 
- p.iR)-pm ' ^''^ 
_ h{R) - MO) .2g. 

^ -p.{R)-P.{oy ^^^^ 

where the reduced kernel of the momentum operator is defined as usually, by P2{R) = 
P2{R)/I{R). Expressions (27) and (28) can be very easily evaluated, especially in view 
of the fact that the momentum kernel can be calculated as a Laplacian of the overlap 
kernel: 

P2{R) = -h^AnHR). (29) 

Of course, within the GOA, the results are exactly the same as those given by the 
zero-point-motion correction and PY mass, Eqs. (10) and (11). However, expressions 
(27) and (28) do not rely on the GOA. They only depend on assuming the quadratic form 
of the Lipkin operator K, Eq. (14). Moreover, variations of exphcit expressions for the 
projected energy -^^(O), hke the ones given by Eqs. (9) or (27), are difficult, while that 
of the Lipkin projected energy (19) can be carried out by the standard self-consistent 
procedure. 

When the quadratic approximation is not sufficient, one can immediately notice 
this fact by a dependence of E^{0) and k on the value of the shift R. In this case, one 
can always switch to higher-order Lipkin operators: 

M 

^ = E hmP'"", (30) 
m=l 

which would require using higher-order Taylor expansions (23), and M different shifts 
Rn, n = 1, . . . , M, instead of one. The only requirement for choosing the shifts -R„ 
is a non-singularity of the matrix A. A dependence of the results on this choice will 
always give a signal that a given order is insufficient. Note that kernels of higher powers 
of the momentum operator can also be calculated in terms of higher derivatives of the 
overlap kernel, in analogy with Eq. (29). However, in view of the fact that for the 
translational symmetry the GOA works so nicely, in this study there does not seem to 
be any immediate necessity to go to higher orders, and the simple quadratic correction 
will suffice. 



2.5. The direct part of ($|p2|$) 

Standard calculations for the SLy4 Skyrme functional [9] are performed for the Lipkin 
operator of K = T/A, which is a fixed factor of A smaller than the one-body kinetic- 
energy operator T. This procedure simply renormalizes the single particle masses as 
m' = m{l — l/A) [15]. Another standard was adopted for the SLy6 and SLy7 functionals, 
whereby the Lipkin operator of Eq. (14) with the exact mass, that is A; = h'^/2mA, was 
used. In Fig. 5, are compared the Lipkin projected energies calculated with the exact 
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Figure 5. (Color online) The Lipkin projected energies (19) calculated in 9 doubly- 
magic spherical nuclei, and plotted relative to the standard energies calculated for the 
SLy4 Skyrme functional [9] . Open and full squares show results obtained for the exact 
[k — fi^ /2mA] and calculated [Eq. (28)] correcting factors, respectively. Solid line 
shows the fit of the volume and surface terms. 
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Figure 6. (Color online) Exact masses (M = mA, open squares) and the PY masses 
calculated after the Lipkin minimization from Eq. (28) (M = /2k, full squares), 
compared with the PY masses of Fig. 4 (full circles) . 
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Figure 7. (Color online) Same as in Fig. 5 but for the chain of lead isotopes around 

208pb. 



and PY, Eq. (28), masses and Fig. 6 shows values of these masses. In all cases, the PY 
masses were calculated by using in Eq. (14) the shift of |-R| = 2fm. 

First one sees that the Lipkin projected energies obtained with the two-body Lipkin 
operators, Eq. (14), differ from those using the standard one-body operator T/A by up 
to 13 MeV. The mass dependence of this difference can be very well described by the sum 
of the volume and surface terms. Therefore, its major part can easily be absorbed in the 
parameters of the Skyrme functional. This confirms conclusions of Ref. [16]. Second, 
one sees in Fig. 6 that the values of the PY masses differ form the exact ones only by a 
few percent. Therefore, the Lipkin projected energies based on these two prescriptions 
for the normalizing constant k differ very little. 

Results shown in Fig. 7 aim at checking whether the differences between the exact- 
mass and PY-mass correcting factors can influence shell effects in masses of lead isotopes. 
One sees, that these differences induce an almost constant shift of the energy, of the 
order of 0.5 MeV. Moreover, the shell effect at ^°^Pb, induced by replacing the Lipkin 
operator kP"^ by T/A, is very small. 

Therefore, based on these results, one might be tempted to consider T/A as a viable 
alternative to kP^. However, in view of the Lipkin symmetry-restoration method, this 
is not the case. This is shown in Fig. 8, where are compared the GOA properties of the 
reduced kernels of and T, 

m^f\i>(R)) 

where the factor in front of the kinetic-energy operator was chosen so that vr is simply 
the kernel of the direct part of the operator P^ [16]. Again, the quality of the GOA is 
here excellent, so the parameters P22 and characterize the kernels very well, while the 
parameters P20 = ($|-P^|$) and ttq = ($|^T|$) give the standard average values. 
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Figure 8. (Color online) Mass-number dependence of the 0-th (lower panel) and 2-nd 
(upper panel) GOA moments of the two-body momentum operator squared and its 
direct part, see Eqs. (31) and (32). 

In the bottom panel of Fig. 8, one sees that the direct parts (ttq) increase almost 
linearly with the particle number A, while the true average values (p2o) increase more 
like A^^^. In a heavy nucleus like ^°^Pb, these trends result in p2o being overestimated 
by TTo by about a factor of 3, which confirms the results obtained in Ref. [17]. Such 
mass dependencies explain why this error can be fairly well absorbed within the volume 
(~y4) and surface (~y4^/^) energies of the Skyrme functional. 

However, in the top panel of Fig. 8, one sees that the second moments of these 
kernels, P22 and 112, behave quite differently. The second moment of the true kernel 
(P22) increases quite fast, like A^^"^, while the direct part (112) only as A^^"^. In ^°^Pb, 
these trends result in P22 being underestimated by 112 by about a factor of 200. It is 
then obvious that one cannot replace the true average values by the direct 

parts, no matter which value of the multiplicative factor is used. 
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One can note here in passing that, within the GOA, Eq. (29) gives the relation 
P22 — fpio' which is very well fulfilled by the numerical results shown in Fig. 8. Than, 
again within GOA, the average values of T in the projected states can be calculated in 
analogy to Eq. (9), and are proportional to 

^GOA/pN ^ 3772 2 (oo\ 

The fact that it2 is a factor of 200 too small in ^osp]^ means that the function 7r(P) is not 
able to flatten the dependence of £'$(P), and thus the method based on replacing 
kP"^ by T/A does not lead to a proper VAP estimate of the projected energy. 

3. Conclusions and perspectives 

In the present study, I have analyzed the translational-symmetry restoration by the 
approximate variation after projection based on the Lipkin method. For translational 
symmetry, the Gaussian Overlap Approximation gives an excellent representation of 
numerical results obtained in doubly-magic spherical nuclei by using the Skyrme energy 
density functionals. The Lipkin method is based on subtracting from the Skyrme energy 
the center-of-mass kinetic energy with the Peierls-Yoccoz mass, and not with the true 
mass. However, calculations show that the Peierls-Yoccoz mass is only a few percent 
different than the true mass of the nucleus. 

I have also studied properties of the direct part of the center-of-mass kinetic 
energy and I showed that this direct part does not fulfill Lipkin conditions and thus 
is not a proper correction within the variation after projection method to restore the 
translational symmetry. In fact, it seems that there is no good argument in support of 
neglecting the exchange part of the center-of-mass kinetic energy at all, apart from the 
fact that it is a cheap method. Nevertheless, shell effects induced by the exchange part 
appear to be weak, and thus the exchange part may only weakly influence the overall 
agreement of the calculated masses with experimental data. 

Due to the fact that the three Cartesian components of the total momentum 
operator commute, the Lipkin method to restore the translational symmetry can be 
easily generalized to deformed nuclei. In this case, one obtains two or three different 
Peierls-Yoccoz translational masses in axially or triaxially deformed nuclei, respectively. 
This natural result reflects the fact that in deformed nuclei the momentum fluctuations 
along the principal axes of the mass distribution are different. 

The methods presented in this study may constitute an interesting alternative to 
the Lipkin-Nogami expressions [7, 18] for calculating the correcting factor A2, required 
for the approximate particle-number-symmetry restoration. Indeed, the Lipkin-Nogami 
expressions are difficult to implement in realistic calculations, and approximate work- 
around procedures have been used in practical approaches [19, 4]. Expressions analogous 
to Eqs. (28) and (29), which are based on the transition energy and overlap kernels with 
respect to a simple shift in the gauge angle, could open here a quite useful new possibility. 
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The largest impact of the Lipkin method can be expected for the rotational- 
symmetry restoration. Here, the Lipkin operator corresponds to the total angular 
momentum squared, which, being included before variation in the functional, may 
induce deformation in all nuclei, including the magic ones. This may create substantial 
rotational corrections in these nuclei, in analogy to those that are obtained by 
minimizing the angular-momentum projected energies after variation over deformation, 
see e.g. Ref. [20]. However, the Lipkin method may magnify these corrections even more, 
because it treats them before variation. Since in the case of the rotational symmetry, 
the Lipkin corrections can be expected to be strongly shell dependent, they may have a 
strong impact on the agreement of the calculated masses with experimental data. 

The same methodology can also be applied to restore the isospin symmetry in 
nuclei. Here, the Lipkin operator corresponds to the total isospin operator squared, 
cf. Refs. [21, 22], and the Peierls-Yoccoz mass must be calculated by considering kernels 
of the nuclear energy only. In this way, the dependence of this nuclear energy on the 
total isospin is flattened, which this leaves to the Coulomb energy the possibihty of 
inducing the correct isospin-mixing effects. 

It is also worth noting that the Lipkin method with the Peierls-Yoccoz collective 
mass is based on properties of the energy kernels for relatively small shifts of the 
collective coordinates. Therefore, singularities of these kernels, see Ref. [23] and 
references cited therein, which plague the exact symmetry-restoration methods within 
the energy-dcnsity-functional approaches, are not causing problems. This observation 
has to be substantiated by analyzing the singularity- free corrected kernels [24, 25, 26] 
and checking that the proposed corrections are appropriately small for small shifts of 
the collective coordinates. 

Altogether, the Lipkin method reviewed in this study can be systematically applied 
to restore, within the mean-field or energy-density-functional theories, all broken 
symmetries. This method constitutes a practical alternative with respect to the exact 
projection techniques, which are very costly and thus cannot be applied to several 
broken symmetries simultaneously. The work on implementing the Lipkin method to 
symmetries other than the translational one studied in this work is now in progress. 

Finally, the Lipkin method can also be used for approximate calculation of average 
values of observables in projected states. However, this method cannot replace the 
real projection in cases when the good quantum numbers are necessary to calculate 
matrix elements and properly account for selection rules. Nevertheless, by projecting 
good quantum numbers from the state that minimizes the Lipkin projected energy one 
probably obtains the best viable alternative to the exact variation after projection. 

This work was supported in part by the Polish Ministry of Science, by the Academy 
of Finland and University of Jyvaskyla within the FIDIPRO programme, and by the 
UNEDF SciDAC CoUaboration under the U.S. Department of Energy grant No. DE- 
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